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ABSTRACT 


“  In  order  to  improve  the  accuracy  of  in-plane  responses  of  shear  deformable  composite  plate 
theories,  a  new  laminated  plate  theory  was  developed  for  arbitrary  laminate  configuration  based  upon 
Reissner's  (1984)  new  mixed  variational  principle.  To  this  end,  across  each  individual  layer,  piecewise 
linear  continuous  displacements  and  quadratic  transverse  shear  stress  distributions  were  assumed.  The 
accuracy  of  the  present  theory  was  examined  by  applying  it  to  the  rylindrical  bending  problem  of  lam¬ 
inated  plates  which  had  been  solved  exactly  by  Pagano  (1969).  A  comparison  with  the  exact  solutions 
obtained  for  symmetric,  antisymmetric  and  arbitrary  laminates  indicates  that  the  present  theory  accu¬ 
rately  estimates  in-plane  responses,  even  for  small  span-to-thickness  ratios. 


1.  INTRODUCTION 


Thick  laminated  plates  and  shells  find  extensive  applications  as  structural  elements.  As  a  result, 
various  approximate  theories  have  been  developed  in  an  effort  to  properly  assess  their  mechanical 
behavior  under  static  and  dynamic  loads. 

In  a  series  of  papers  Paganc  (1969,  1970,  1971,  1972)  gave  the  exact  solution  for  the  problem  of 
cylindrical  bending  and  simply  suppo  ted  plates.  He  pointed  out  the  importance  of  considering  shear 
deformation  effects  m  order  to  accurately  predict  the  plate  lateral  deflection  and  the  necessity  of 
improving  the  thickness  variation  of  in-plane  displacements,  which  are  assumed  to  be  C’-linear  func¬ 
tions  in  both  classical  plate  theory  (CPT)  and  Reissner-Mindlin  plate  theory  (FSD). 

One  of  the  earliest  attempts  in  deriving  an  approximate  theory  for  laminates  is  credited  to  Yu 
(19591  He  investigated  the  plane  strain  problem  of  isotropic  sandwich  plates  by  assuming  piecewise 
linear  displacement  distributions  Yang.  Norris  and  Stavsky  (1966)  extended  Reissner-Mindlin  plate 
theory  to  the  case  of  an  arbitrar>  number  of  bonded  anisotropic  layers.  Whitney  and  Pagano  0970) 
using  this  latter  approach,  and  later  Whitney  (1972).  concluded  that  the  introduction  of  shear  deforma¬ 
tions  cannot  improve  the  in-plane  stress  distributions  as  determined  from  classical  plate  theory. 

As  a  remedy  to  these  difficulties,  higher  order  theories  have  been  proposed  in  which  the  displace¬ 
ment  assumptions  are  expressed  in  terms  of  power  series  in  the  thickness  variable.  The  number  of 
plate  equations  of  such  theories  does  not  increase  with  the  number  of  layers  Theories  including  qua¬ 
dratic  variations  (V^liitney  and  Sun,  1973.  and  Nelson  and  Lorch.  1977)  and  cubic  variations  (Hilde¬ 
brand.  Reissner  and  Thomas.  1949.  and  Lo,  Christensen  and  Wu.  1977)  of  in-plane  displacements 
through  the  plate  thickness  belong  to  this  category.  Reddy's  (19841  high-order  theory  is  obtained  by 
imposing  the  condition  of  vanishing  transverse  shear  strains  on  the  top  and  bottom  surfaces  of  the 
plate.  N^Tiitney  (1972)  derived  in-plane  displacements  by  integrating  the  transverse  shear  strains 
deduced  by  Whitney  and  Pagano  (1970).  This  resulted  in  a  higher  order  approximation  which  accu¬ 
rately  predicted  in-plane  strains,  but  the  resulting  modified  stresses  did  not  necessarily  satisfy  the  origi¬ 


nal  plate  equilibrium  equations 
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Another  class  of  higher  order  theories  employs  the  displacements  to  be  continuous  and  piecewise 
smooth  functions,  i  e..  smooth  within  each  layer  (Yu,  1959)  Durocher  and  Solecki  have  followed  this 
approach  to  analyze  transversely  isotropic  sandwich  plates.  Mau  (1973),  Srinivas  (1973)  and  Seide 
(1980)  have  considered  the  case  of  an  arbitrary  number  of  layers.  No  shear  correction  factors  were 
inir;)duced  by  the  last  two  authors.  In  all  these  theories,  the  number  of  field  equations  and  edge  boun- 
dar\  conditions  depends  on  the  number  of  layers. 

Most  of  the  theories  discussed  so  far.  which  are  classified  as  displacement-based  theories,  suffer 
from  a  common  deficiency:  constitutive  equations  lead  to  discontinuous  interlaminar  stresses.  This 
shortcoming  has  been  overcome  in  the  theory  proposed  by  Murakami  (1985)  which  is  based  upon 
Reissner's  (1984)  new  mixed  variational  principle  This  theory  is  obtained  by  superposing  a  zig-zag 
shaped  -linear  function  to  the  in-plane  displacements  given  by  Reissner-Mindlin  plate  theory.  It  was 
later  extended  by  including  Legendre  polynomials  in  the  displacement  variations  across  the  plate  thick¬ 
ness  (Toledano  and  Murakami.  1985)  In  both  these  theories,  the  number  of  equilibrium  equations 
and  edge  boundary  conditions  are  independent  of  the  number  of  layers.  Thus,  the  theories  have  some 
limitations  when  applied  to  composite  plates  with  arbitrary  laminate  configuration. 

In  order  to  improve  the  in-plane  response  of  composite  plates  for  arbitrary  laminate  configuration. 
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a  new  theory  which  accounts  for  transverse  shear  deformations,  has  been  developed  by  assuming  piece- 
wise  continuous  in-plane  displacement  distributions.  In  order  to  guarantee  continuity  of  interlaminar 


stresses.  Reissner's  (1984)  new  mixed  variational  principle  has  been  invoked  by  taking  the  transverse 
stresses  to  be  quadratic  functions  of  a  local  thickness  coordinate  across  each  layer.  Governing  equa¬ 
tions  and  consistent  boundary  conditions  are  then  deduced.  The  advantage  of  using  Reissner's  new 
mixed  variational  principle  is  that  it  automatically  yields  the  appropriate  shear  correction  factors  for  the 


transverse  shear  constitutive  equations  A  comparison  with  Pagano's  (1969)  exact  elasticity  solution  for 
symmetric,  antisymmetric  and  arbitrary  laminates  in  cylindrical  bending,  shows  that  the  proposed  theory 
can  accurately  predict  in-plane  displacements  and  stresses  at  low  span-to-thickness  ratios.  To  further 
assess  ihe  range  of  applicability  of  the  present  theory,  results  previously  obtained  by  Murakami  (1985) 


are  also  presented 


2.  FORMILATION 


Consider  an  N-layer  laminated  composite  plate  of  uniform  thickness  h,  as  shown  in  Fig.  1.  A 
cartesian  coordinate  system  is  chosen  such  that  the  middle  surface  of  the  plate  occupies  a  domain  D  in 

the  .V],  x;-plane,  the  xj-axis  being  normal  to  this  plane.  The  following  notation:  (  )•*’,  k  »  1,2 . .V 

will  designate  quantities  associated  with  the  /r'*-layer.  The  thickness  of  each  layer  is  in  which 

the  \olume  fractions  n'*'  satisfs  the  relation 


I  n'*' 
c-i 
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(1) 


Unless  otherwise  specified,  the  usual  cartesian  indicial  notation  is  employed  where  latin  and  greek 
indices  range  from  1  to  3  and  1  to  2.  respectively  Repeated  indices  imply  the  summation  convention 
and  I  )  is  used  to  denote  partial  differentiation  with  respect  tox, . 

With  the  help  of  the  foregoing  notation,  the  governing  equations  for  the  displacement  vector  u/*' 
and  stress  tensor  o-,',* '  associated  with  the  A"' -layer  are: 

a)  Equilibrium  Equations 

<>.^“'-0  ;  (2) 

where/  are  body  forces. 

b)  Constitutive  Equations  for  Monoclinic  Layers 
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where  C,.  are  the  elastic  constants  and  C„  iij  -  1.2.6)  represent  the  reduced  stiffnesses  introduced  b\ 


Whitney  and  Pagano  (1970); 


c)  Strain-Displacement  Relations 


d)  Interface  Continuity  Conditions 

^  u  ^  u-ii  ^  u  I  _  ^u.i.  .  f.  _  j  2,  .  •  ■  ,  AT-i 

e)  Upper  and  Lower  Surface  Stress  Conditions 

<73',''-  T;*  on  ^3”  7 

(T3?'-rr  on  X3---| 
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(5) 


(6a) 

(6b) 


The  objectives  in  developing  the  present  composite  plate  theory  are  :  first,  to  improve  the 
assumed  variations  of  in-plane  displacements  through  the  plate  thickness;  second,  to  take  into  account 
the  effect  of  transverse  shear  deformation  and  derive  appropriate  constitutive  equations;  and  finally,  for 
arbitrary  laminate  configurations,  to  assess  and  compare  its  range  of  applicability  with  the  laminated 
plate  theory  proposed  by  Murakami  (1985). 

In  order  to  carry  out  this  task,  Reissner’s  (1984)  new  mixed  variational  principle  for  displace¬ 
ments  and  transverse  stresses  was  applied  to  the  .V -layer  composite  plate. 
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where  dDr  denotes  the  boundary  of  D  with  outward  normal  on  which  tractions  T,  are  prescribed 
and  /4'‘'  is  the  jr;-domain  occupied  by  the  A"" -layer.  .Mso  rJ,*'  denote  the  approximate  transverse 


3.  TRIAL  DISPLACEMENT  AND  TRANSVERSE  STRESS  FIELDS 


The  present  laminated  plate  theory  which  accounts  for  transverse  shear  effects  is  obtained  by 
assuming  a  linear  variation  of  the  in-plane  displacements  across  each  individual  layer,  as  shown  in 
Fig.  2.  Transverse  displacements  are  taken  to  be  constant  throughout  the  entire  thickness  of  the  plate. 
Therefore  shear  strains  are  constant  within  each  layer,  but  differ  from  layer  to  layer.  This  approach  has 
been  previously  adopted  by  Seide  (1980)  for  the  .V-layer  case.  However  no  mention  of  shear  correc¬ 
tion  factors  was  made  in  his  paper. 

The  appropriate  trial  functions  used  in  connection  with  Reissner's  mixed  variational  principle 
Eq  (7)  are  chosen  to  be; 

a)  Trial  Displacement  Field 

'  (.V, )  =  L'a*  ■ ' '  (.Vg  ‘  (.V 3* ' )  -i-  '  (.x-^  )gi‘'  ’  (.vi* ' )  (8a) 

U3*'(.x,)-  fjCx,,)  (8b) 

where  gj* *  (x j* ' )  =  4  (- 1  )°* ‘  -^rr:  (9) 

2  h 

xj' '  IS  a  local  x^-coordinate  system  with  its  origin  at  the  center  xjo'  of  the  -layer,  i.e. 

xj*’  =  xj- .vjg'  (10) 

From  Eqs.  (8a)  and  (9)  it  is  seen  that  La‘’(A  -  1.2,  ■  ..V-D.  and  Lj''  represent  the  values 

of  at  the  interface,  top  and  bottom  surfaces  of  the  plate,  respectively  Also.  Eq.  (8)  satisfy  inter¬ 
face  displacement  continuity  conditions  Eq  (5a). 

h)  Trial  Transverse  Stress  Field 


where 


■j^'(x  )-  0:‘’(x,)F,(.-)^  C-"(x,tf3(.-)^  r'‘’(x3)F3(--) 
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n'”  h 

Also  0i"=/Tli'dx3  (14) 

>4“' 

In  Eq  (11a).  Tj*'"  and  are  the  values  of  '  at  the  top  and  bottom  surfaces  of  the  k"'  layer, 
respectively.  From  Eq.  (6)  one  has 

and  r/'’=r,-  '15) 

Eq  (in  satisfy  the  interface  stress  continuity  conditions  Eq  (5b)  Due  to  the  approximation  for  nj*' 
which  yields  eiV  =  0,  0-33’  becomes  a  reactive  stress.  Consequently,  o-j^'  can  be  determined  by 


integrating  the  third  equilibrium  equation. 
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4.  LAMINATED  PLATE  EQUATIONS 

Substituting  Eqs  (8)  and  (11)  into  Eq.  (7)  and  using  Gauss'  theorem,  one  obtains; 
a)  Equilibrium  Equations 
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b)  Boundary  Conditions 
Specify  '  or 
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Eq  (18)  constitute  (2.V  +  3)  conditions.  It  can  be  seen  that  the  natural  edge  traction  boundary  condi¬ 
tions  are  coupled,  i.e.  force  and  moment  resultants  acting  on  two  adjacent  layers  are  involved, 
c)  Constitutive  Equations  for  Monoclinic  Layers 
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In  Eqs  (19a. b)  k  ranges  from  1  to  N,  while  in  Eqs  (19c. d)  k  ranges  from  1  to  (,V-1).  Eqs.  (19)  can 
be  solved  for  ’  and  Tj* '  in  terms  of  L'r  ’  and  As  a  result,  the  quantities  can  be  deter¬ 
mined  as  functions  of  these  displacement  variables.  Such  expressions  will  automatically  include 
appropriate  shear  correction  factors  by  virtue  of  Reissner's  mixed  variational  principle. 

The  remaining  constitutive  equations  for  ,Vq^’  and  are  obtained  b>  substituting  (3a).  (4) 
and  (8)  into  ( 1 7a)  to  vield 
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5.  CYLINDRICAL  BENDING  OF  LAMINATED  PLATES 

In  order  to  lest  the  accuracy  and  assess  the  range  of  applicability  of  the  present  theory,  cylindrical 
bending  of  composite  plates  under  sinusoidal  loading  is  considered.  The  plate  is  simply  supported  at 
the  ends  .V|  =  0  and  /,  and  is  infinitely  long  in  the  jr^-direction.  The  prescribed  boundary  conditions  on 
the  top  and  bottom  surfaces  of  the  plate  are 

r*  =  0  ,  Tf  =  <7  sin  on  .V3  =  (21a) 


rr  -  rj  =  0 


on  X)  = 


h 

2 


(21b) 


The  boundary  conditions  for  the  simply  supported  ends  are.  from  (18) 

L'3=  0  at  .V,  =  0./ 


(.vir  +  ‘ 

X'  U-  111 

1 

1 

Nil  1 

/r 

„ixi  „u-i.  1 

=  0  at  .v,  =  0.1 


(22a) 


(22b) 


For  simplicit\.  only  cylindrical  bending  of  laminated  plates  consisting  of  orihoiropic  layers  will  be  con¬ 
sidered  in  the  subsequent  analysis.  However,  no  additional  restrictions  are  imposed  on  layer 
thicknesses,  elastic  moduli  and  stacking  sequence.  In  this  case,  there  holds 


C),'* '  »  ’  =  0  (23) 

C\V  =  0  and  -  1  (24) 

.^5  a  result,  terms  involving  and  will  drop  out  from  Eqs  (19a. c'.  while  Eqs  (19b.d)  will  not 

appear  altogether  The  remaining  equations  for  Ql‘''  and  Fi*  can  be  written  in  matrix  form  as 


Qt- p 

(25a) 

/iif.ir,  -  0 

(25b) 

re 


and  [fill.  [Dll  and  [fil  are  matrices  of  dimensions  .V  x  (,V-1)  x  S  and  (A'-l)  x  (V-l). 

respectively.  These  matrices  are  only  functions  of  n'"’  and  The  righi-hand  side  of  Eq.  (25a) 

contains  the  displacement  variables  L'/*’  and  Eqs.  (25)  can  be  solved  for  Ql‘‘'  and  T]'*'  to  yield 


/I  r,  =  -  [Tr'l[D,l^i  (27a) 

^1-  ([/!+  [fiil[fr'][D,l)-' *  (27b) 

where  [/I  is  the  .V  x  .V  identity  matrix. 

The  remaining  constitutive  equations  are  simply,  from  (20) 

.Vi'f  =  (28a) 

•V/f '  -  '  jt’iV"  -  i’iV)  (28b) 

,V<V  =  .v'V  =  =  0  (28c) 


Using  surface  boundary  conditions  (21).  the  equilibrium  equations  (16)  for  cylindrical  bending  in  the 
absence  of  body  forces  reduce  to 
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(29b) 
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■'31.1 


q  sin 


I 


(29d) 


The  form  of  the  dependence  on  the  displacement  variables  U;*'  and  L  y  of  the  constitutive  equations 
(27b)  and  (28a. b).  and  the  nature  of  the  applied  load  suggests  the  folloivtng  expressions  for  the 


displacements 


t/,'*’  -  h  L','* ’  cos  and  h  Uy  sin  (30) 

where  and  Vy  are  nondimensiona)  quantities  by  definition.  It  is  easily  proven  that  the  boundary 
conditions  (22'  are  satisfied  when  (30'  are  substituied  therein 

Finallv,  inserting  (30)  into  the  coistitutive  equations  <27b)  and  (28a.b).  and  these  in  turn  into 
the  equilibrium  equations  (29),  yields  a  system  of  (S+2)  algebraic  equations  with  the  (A'  +  2)  nondi- 
mensional  amplitudes  t','*'  and  L'y  as  unknowns.  This  system  is  conveniently  written  in  matrix  form 
as 

[X]  0~  F  (31) 

where  .  •  •  •  .  L',''’  ,  (32a) 

T=[0.0 . O.q]’’  (32b) 

and  (A’]  is  a  (,V  +  2)  x  (A'  +  2)  matrix. 


6.  NUMERICAL  RESULTS  AND  DISCUSSION 


In  order  to  test  the  accuracy  of  the  present  theory,  the  problem  of  cylindrical  bending  of  an 
infinitely  long  strip  under  sinusoidal  loading  is  re-examined.  The  exact  elasticity  solution  has  been 
given  by  Pagano  (1969),  where  a  symmetric  three  layer  cross-ply  laminate  was  considered,  the  0°  layers 
being  at  the  outer  surfaces  of  the  plate.  The  elastic  properties  are 

c  c 

for  the  0°  layers  -  25.062657  .  -  0.5  (33a) 

Ej  Ef 

C  C 

and  for  the  90°  layers  *  1.002506  .  -=^-0.2  (33b) 

Lt  Ej 

where  Ej  is  a  reference  modulus. 

Following  Pagano's  (1969)  nondimensionalization.  the  displacements  and  stresses  are  written  in 
the  form 

u,  -  ^  /,  •  -  Y  ^2 

(34) 

o^ii  -  ^  <^11  ’■j  • 


•’^3  A  C  I 

—  and  S  =  T 
h  h 


In  the  various  curves  the  solid  line  represents  the  exact  solution  while  the  results  of  the  present 
theory  are  shown  by  a  broken  line.  Also  shown,  for  comparison  purposes,  are  the  results  given  by  the 
laminated  plate  theory  proposed  by  Murakami  (1985)  which  are  represented  by  a  dashed-dotted  line.  It 
is  a  shear  deformable  theory  obtained  by  superposing  to  the  linear  variations  of  the  Reissner-Mindlin 
theory  a  zig-zag  in-plane  displacement  variation  across  the  plate  thickness.  For  brevity .  this  theory  will 
be  called  here  "The  First-Order  Zig-Zag  Theory*'  and  abbreviated  as  ZZ 

For  a  symmetric  3-layer  cross-ply  laminate  (0/90/0)  with  layers  of  equal  thickness.  Table  1  shows 
the  values  of  the  central  deflection  uj  obtained  from  the  different  theories  for  a  span-to-thicknes<>  ratio 
S  -  4  As  observed,  the  present  theory  and  the  First-Order  Zig-Zag  theory  yield  exactly  the  same 


and 


numerical  result  This  is  also  true  for  the  thickness  variations  of  the  in-plane  displacement  Cj*' 
normal  stress  As  shown  in  Figs  3a. b,  very  close  agreement  is  found  between  both  theories  and 
Pagano  s  exact  solution  Jt  should  be  pointed  out  that  the  numerical  results  obtained  by  Seide  (1980) 
differ  slightly  from  those  given  here  since  no  shear  correction  factors  were  introduced  by  him 

The  present  theory  was  next  tested  for  a  symmetric  5-layer  cross-ply  laminate  (0/90/0/90/0)  with 
layers  of  equal  thickness.  The  central  deflection  uy  for  S  -  4  is  shown  m  Table  1.  where  closer  agree¬ 
ment  between  the  present  theory  and  the  exact  solution  is  obserxed  as  compared  to  the  first-order  zig¬ 
zag  theory.  The  distributions  across  the  plate  thickness  of  in-plane  variables  ui'*'  and  an  '  are  com¬ 
pared  in  Figs  4a.b.  It  is  seen  that  the  present  theory  has  improved  upon  ZZ.  especially  in  the  interior 
layers  of  the  plate. 

.^n  antisymmetric  4-layer  cross-ply  laminate  (0/90/0/90)  with  layers  of  equal  thickness,  was  also 
examined  In  this  case,  the  error  in  the  central  deflection  between  the  exact  solution  and  ZZ  is  quite 
large  (-•  2T''o).  while  the  present  theorx  still  gives  a  satisfactory  value  (see  Table  1).  The  variations 
across  the  plate  thickness  of  the  in-plane  displacement  «|“'  and  normal  stress  Jrn'  are  shown  in 
Figs  5a. b  for  5  =  4.  respectively.  From  the  curves  for  Cj*'.  it  is  seen  that  the  first-order  zig-zag 
theory  deviates  significantly  from  the  exact  solution  at  the  bottom  layer  of  the  plate.  On  the  other 
hand,  the  present  theory  is  in  very  good  agreement  with  the  exact  solution. 

To  further  assess  the  range  of  applicability  of  both  the  present  theory  and  the  first-order  zig-zag 
model,  arbitrary  laminate  configurations  consisting  of  3.  4  and  5  layers  were  tested.  Three  different 
materials  were  used,  with  the  following  elastic  properties 

c  c 

Material]  :  -=^=  1.002506  .  -r^  -  0.2  (36a) 

tj 

Material  2  :  -  32.631  .  -  8.21  (36M 

Et  Et 

Material  3  :  -r^  -  25  062657  .  =  o  5 

Et  Et 


( 36c ' 
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The  laminate  configurations  corresponding  to  the  three  cases  examined  are  shown  in  Table  2.  In  all 
cases.  S  “  4. 

The  values  of  the  central  deflection  for  the  three  laminate  configurations  are  given  in  Table  3.  It 
can  be  observed  that  the  discrepancies  with  the  exact  solution  are  larger  in  the  case  of  the  first-order 
zig-zag  theory  than  in  the  case  of  the  present  theory.  In  particular,  the  error  in  Uj  for  A’  —  5  is  as  high 
as  49'’'(.  The  variations  across  the  plate  thickness  of  the  in-plane  displacement  ul*'  and  normal  stress 
5)1'  are  shown  in  Figs  6.  7  and  8  for  .V  «  3.4  and  5.  respectively  As  expected,  the  symmetric  distri¬ 
butions  obtained  by  the  present  theory  and  ZZ  for  .V  -  3  and  5  no  longer  hold  for  arbitrary  lamnate 
configurations.  In  all  the  cases  considered,  the  present  theory  is  still  in  good  agreement  with  the  exact 
solution,  except  possibly  at  the  top  layer  of  the  4-ply  laminate  (see  Fig  7).  On  the  other  hand,  it  is 
seen  that  the  first-order  zig-zag  theory  deviates  considerably  from  the  exact  solution. 

Some  important  points  are  now  discussed  The  accuracy  of  the  present  theory  can  be  improved  by 
dividing  each  ply  into  a  finite  number  of  sub-layers,  but  at  the  expense  of  increasing  computer  storage 
Also,  the  proposed  theory  possesses  two  main  drawbacks,  which  are  .  first,  the  number  of  equilibrium 
equations  and  edge  boundary  conditions  increases  with  the  number  of  layers;  and  second,  due  to  the 
coupling  of  the  natural  edge  boundary  conditions  (see  Eq  (18)).  no  clear  physical  meaning  seems  to  be 
associated  with  them. 

Next,  as  a  possible  explanation  for  the  increased  discrepancies  in  the  case  of  arbitrary  laminates 
between  the  exact  solution  and  the  first-order  zig-zag  theory,  the  following  argument  is  proposed  The 
inclusion  of  the  zig-zag  shaped  C^-function  was  motisated  by  the  displacemeni  microstruciue  of 
periodic  laminated  composites  (Murakami.  Maewal  and  Hegemier.  1981)  Obviously,  for  arbitrary  lam¬ 
inate  configurations,  this  periodicity  is  destroyed.  Therefore,  the  first-order  zig-zag  theory  should  be 
expected  to  break  down  in  these  particular  cases. 

Finally,  it  is  worth  mentioning  that  an  advantage  of  the  zig-zag  model  is  that  the  number  of  equa¬ 
tions  to  be  solved  for  the  displacement-type  variables  are  independent  of  the  number  of  layers  This  is 
certainly  not  true  for  the  present  theory.  In  design  problems  incoUmg  laminates  with  a  large  number 
of  layers,  these  considerations  should  be  taken  into  account 
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7.  CONCLUSION 

A  composite  plate  theory,  which  accurately  predicts  in-plane  responses  of  arbitrary  laminates,  was 
developed  by  assuming  a  linear  variation  of  in-plane  displacements  and  a  quadratic  variation  of 
transverse  stresses  across  each  individual  lamina.  Transverse  displacements  were  kept  constant 
throughout  the  entire  plate  thickness  Governing  equations  and  appropriate  boundary  conditions  were 
then  deduced  from  Reissner's  (1984)  new  mixed  variational  principle  The  accuracy  of  the  theory  was 
examined  for  the  case  of  cylindrical  bending  of  an  infinitely  long  strip  and  compared  with  the  exact 
solution  given  by  Pagano  (19691.  The  results  obtained  by  the  laminated  plate  theory  proposed  by- 
Murakami  (1985).  and  called  "The  First-Order  Zig-Zag  Theory"  were  also  showTi.  Values  for  the  central 
deflection  and  in-plane  displacement  and  normal  stress  for  symmetric,  antisymmetric  and  arbitrary  lam¬ 
inates  were  presented  It  was  observed  that  the  first-order  zig-zag  theory  gave  satisfactory  results  only 
for  symmetric  cross-ply  laminates.  The  discrepancies  between  this  theory  and  the  exact  solution  were 
more  pronounced  for  laminates  of  arbitrary  configuations  than  for  antisymmetric  cross-ply  laminates. 
On  the  other  hand,  in  all  the  cases  considered,  the  present  theory  was  in  close  agreement  with  the  exact 
solution  Despite  certain  shortcomings,  the  good  correlation  with  the  exact  elasticitv  solution  indicates 
that  the  proposed  theory  may  prove  useful  in  the  investigation  of  the  mechanical  response  of  arbitrary 
laminates. 
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Table  1.  Centra)  Deflection  Uj  for  Symmetric  3  and  5-Layer  and  Antisymmetric  4-Layer  Cross- 
Ply  Laminates  (S  »  4) 


(0/90/0) 

5 

(0/90/0/90/0) 

/V-  4 

(0/90/0/90) 

Exact  Solution 

2.887 

3.044 

4.181 

Present  Theory 

2.907 

3.059 

4.202 

l®*-Order  Zig-Zag 

2.907 

3.018 

3.316 

-  2  - 

Table  2.  Geometric  and  Physical  Properties  of  Arbitrary  3.  4  and  5-Layer  Laminates  (S  -  4) 
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Table  3.  Central  Deflection  U3  for  Arbitrary  3,  4  and  5-Layer  Laminates  (S  -  4) 


3 

A'  -  4 

A-  5 

Exact  Solution 

2.341 

1.665 

2.456 

Present  Theory 

2.364 

1.644 

2.467 

l^'-Order  Zig-Zag 

1.992 

1.303 

1.261 
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